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Brillouin scattering has been used to measure the dynamic structure factor of glycerol as a function 
of temperature from the high temperature liquid to the glassy state. Our investigation aims at 
understanding the number and the nature of the relaxation processes active in this prototype glass 
forming system in the high frequency region. The associated character of glycerol is reflected by a 
rather simple relaxations pattern, while the contributions coming from intra-molecular channels are 
negligible in the GHz frequency region. The temperature behavior of the characteristic frequency 
and lifetime of the longitudinal acoustic modes is analyzed, suggesting that a phenomenological 
model which only includes the structural (a) process and the unrelaxed viscosity is able to catch the 
leading contributions to the dynamics of the density fluctuations. This ansatz is also supported by a 
combined analysis of light and inelastic x-ray scattering spectra. The temperature dependence of the 
characteristic time of the a-process, r a , obtained by a full-spectrum analysis conforms to the a-scale 
universality, i.e. the values r a revealed by different experimental techniques are proportional the ones 
to the others. The non-erogodicity parameter smoothly decreases on increasing the temperature, 
and no signature of the cusp-like behavior predicted by the idealized mode coupling theory and 
observed in other glass-formers is found in glycerol. 



I. INTRODUCTION 

Density fluctuations hold a key role in the dynam- 
ics of glass-forming systems: the glass formation itself 
is related to the progressive arrest of long time density 
fluctuations. The characteristic time t q of these fluc- 
tuations exhibits a steep increase for temperatures ap- 
proaching the glass transition temperature, T g , ranging 
from ^-T0~ 12 s well above the melting point to ~10 2 s 
at T g . For this reason more than one technique is usu- 
ally required to cover the whole time or frequency range 
where this dynamics takes place. Inelastic neutron, x-ray 
and light scattering are the experimental techniques most 
widely used to study the density fluctuations in super- 
cooled systems. Among these, Brillouin light scattering 
(BLS) is a powerful tool to gain access to the spectrum 
of density fluctuations in the GHz region. 

In the earliest BLS studies of glass forming systems, 
the data analysis was focused on measuring the fre- 
quency shift and linewidth of the Brillouin peaks, i.e. the 
phase velocity and absorption of the longitudinal acoustic 
modes 0, 0, yj . The analysis of the temperature depen- 
dence of these parameters was used to determine either 
the temperature dependence of the relaxation time (ac- 
tivation plot) once the shape of the relaxation function 
was fixed, or, conversely, the shape of the structural re- 
laxation (master plot) once the r(T) behavior was given 

!. 

More recently, the improvements of interferometric se- 
tups allowed the experimentalists to obtain high-contrast 
and high-resolution BLS spectra in a wider frequency 
window ranging from few hundreds of MHz to few tens 
of GHz, thus gaining information on the shape of both 
the Brillouin peaks and the quasi-elastic central region 



known as Mountain peak 5J. Thanks to these advances, 
it has become possible to carry out full spectrum analy- 
ses in a window of l-=-2 decades in frequency, using ex- 
pressions for the spectrum of density fluctuations derived 
within the Generalized Hydrodynamics formalism Q or, 
equivalently, from a generalized Langevin equation 0. 
The spectrum is thus written in terms of static correla- 
tors and of a frequency-dependent longitudinal modulus 
or, equivalently, longitudinal viscosity. In the earliest 
investigations of supercooled systems, reasonable fits of 
the measured spectra were obtained considering the a- 
relaxation as the only relevant relaxational process. In 
those investigations, the a-relaxation was typically de- 
scribed by the Cole-Davidson (CD) expression, like in 
the case of salol ||, PC @, ZnCl 2 GH EH The main 
limit of this approach is to ignore the possible presence 
of additional processes in the BLS frequency range, faster 
than the CD relaxation. 

A contribution to the fast dynamics of density fluctu- 
ations in supercooled systems is predicted by the mode 
coupling theory (MCT) [12J as the fast component of the 
structural relaxation, also known as the mode coupling (3 
region. During the last decade an important contribution 
to the comparison of the MCT predictions to Brillouin 
spectra has been given by Cummins and coworkers |13| . 
Another important contribution to the fast dynamics of 
density fluctuations in molecular liquids comes from in- 
ternal thermal relaxations These are connected to 
the exchange of energy between acoustic waves and in- 
tramolecular degrees of freedom, like vibrations and ro- 
tations of groups of atoms. It is interesting to notice 
that the Generalized Hydrodynamic approach to Bril- 
louin scattering was originally introduced in the sixties 
by Mountain and Zwanzig p|, to explain exactly the 
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thermal relaxation in Kneser-type liquids |LJ], although 
its importance has been largely underestimated in more 
recent studies of supercooled systems. A clear signa- 
ture of the existence of this relaxation channel in BLS 
spectra can be found in the almost Arrhenius tempera- 
ture dependence of the relaxation time obtained by the 
single-relaxation analysis of different glass forming sys- 
tems, both simple [9| and polymeric (see for instance, 
HI, llll) in nature. 

The main reason why a two-relaxation (a+thcrmal) 
model was not employed for long time in the analysis 
of BLS spectra can be found in the narrowness of the 
frequency window available to this technique. A step 
forward in the solution of this problem has been recently 
proposed, based on the combined use of ultrasonic, BLS 
and inelastic x-ray scattering (IXS) spectra. IXS is a 
novel technique which directly probes density fluctua- 
tions in the high q limit of some nm _1 , extending the 
frequency range accessible to Brillouin scattering inves- 
tigations to the THz region. It is then now possible to 
carry out a full spectrum analysis in a reasonably wide 
spectral region, including both structural, thermal and 
unrelaxed viscosity ("instantaneous") processes. Using 
this method in polybutadiene it was possible to sep- 
arate the contribution of the a relaxation from that of the 
fast intramolecular process, and to reconcile the positive 
tests of MCT previously obtained by neutron scattering 
with BLS results. In ortho-terphenyl, the presence of a 
vibrational relaxation has been confirmed by BLS 

measurements on the single crystal |2l) , and by molecular 
dynamics simulations [22|. 

In this work, we study the high frequency dynamics 
of the glass former glycerol, an associated liquid which 
forms a H-bond network, making it intermediate between 
fragile and strong systems. Different from Kneser liq- 
uids, in associated liquids the ratio between the acous- 
tic absorption a and the classical absorption a c i (only 
due to shear viscosity and thermal conductivity) is very 
small (typically less than 3) and temperature indepen- 
dent. This feature is traditionally attributed to a negli- 
gible coupling of the acoustic waves with internal degrees 
of freedom. This simple relaxation pattern suggests glyc- 
erol as a natural candidate for testing glass formation 
theories. 

Glycerol has already been the object of several stud- 
ies concerning its acoustic properties and dynamical pro- 
cesses. Ultrasonic (US) measurements in the MHz fre- 
quency region have been reported of the adiabatic speed 
and of the absorption of sound H IM HI El • Photon 
correlation spectroscopy (PCS) has been used to charac- 
terize the structural relaxation close to T g [26l |27| . Im- 
pulsive stimulated Brillouin and thermal light scattering 
(ISBS and ISTS) experiments [H H<| have provided in- 
formation on the sound dispersion and absorption and 
on the structural relaxation at low temperatures and for 
different q values. Stimulated Brillouin gain (SBG) spec- 
troscopy has been also used to probe the dynamics of 
supercooled and glassy glycerol 30]. Very recently, the 



density response of supercooled glycerol to an impulsive 
stimulated thermal grating has been studied in the 200- 
340 K temperature range in order to separate the struc- 
tural relaxation from thermal diffusion processes |3l| . 
Rayleigh-Brillouin spectra at very high resolution have 
been presented |32| to reveal acoustic anomalies in glassy 
glycerol at low temperatures. The temperature behavior 
of the structural relaxation process has also been studied 
by depolarized light and neutron scattering (3j| |3 HH i 
and dielectric spectrosco py l3 r^ . l37l l3S| . Moreover, inelas- 
tic x-ray scattering |39l I40 u4lj has been used to inves- 
tigate the high-frequency dynamics of supercooled glyc- 
erol and the nature of the acoustic attenuation at low 
temperatures. Despite of this extensive literature, there 
is no light scattering investigation which covers the en- 
tire range comprised between the high temperature liquid 
and the glassy state. 

Here we report high resolution Brillouin spectra cov- 
ering nearly two decades in frequency in a wide temper- 
ature range from 47 to 441 K. The main purpose of this 
paper is to contribute to the identification of the relax- 
ation processes active in the GHz frequency region of 
glycerol, and to the characterization of their tempera- 
ture evolution. In order to constrain the hydrodynamic 
model proposed for the density fluctuations, we perform 
a combined light and inelastic x-ray analysis. Moreover, 
a comparison is also reported with the results obtained by 
different techniques, including ISTS, ISBS, SBG, dielec- 
tric spectroscopy, photon correlation spectroscopy, and 
depolarized light scattering, in order to test the degree of 
universality of the temperature dependence of the char- 
acteristic time and of the shape parameters of the a re- 
laxation. 

Finally, a preliminary comparison with the predictions 
of the mode coup ling theory of supercooled liquids is 
presented. MCT [TJ, considered the first microscopic 
model of the glass transition, describes the slowing down 
of the structural dynamics in terms of a non-linear cou- 
pling between density fluctuations which, if we neglect 
hopping processes, causes the structural arrest of the sys- 
tem at the critical temperature, T c . Several studies have 
been previously performed to test the MCT predictions 
on glycerol. Depolarized li ght scattering [33, yj, [35| and 
neutron [33 and dielectric |36l I3H I38L l42| data have been 
subjected to MCT checks. All together these studies pro- 
vide a variegated picture of the existence and the eventual 
location of the critical temperature, T c : the different T c 
evaluations span a range comprised between 250 and 310 

k m0msni[3i,[i2. 

One of the main predictions of the idealized mode cou- 
pling theory concerns the temperature behavior of the 
non-ergodicity factor f q — the normalized amplitude of 
the structural relaxation — which should exhibit a cusp 
at T c . In this paper, by measuring the limiting low and 
high frequency adiabatic sound velocity in a wide tem- 
perature range, we derive the non-ergodicity factor |43| 
for the density fluctuations, and look for the possible ex- 
istence of the crossover temperature predicted by MCT. 
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II. THEORETICAL BACKGROUND 



A. Dynamic structure factor 



modulus (or viscosity). In particular, the complex fre- 
quency dependent modulus M(ui) — M'(lo) — iM"(ui) 
can be written as: 



An expression for the dynamic structure factor, and 
thus for isotropic Brillouin spectrum which is propor- 
tional to it, can be derived from the equation of motion 
of the microscopic density p(r, t) . In the limit of low ther- 
mal conductivity, density fluctuations are decoupled from 
temperature fluctuations, and the equation of motion of 
the g-th component of the microscopic density, p(q,t), is 
derived within simple hydrodynamics and comes out to 
be that of a damped harmonic oscillator (DHO): 



q 2 J dt 2 VL dt 



M 



p{q,t) = Q 



(1) 



where pm is the average mass density, r)L the longitu- 
dinal viscosity and M the longitudinal elastic modulus. 
The dynamic structure factor, S(q, u>), obtained from the 
Laplace-Fourier transform of Eq. (JIJ, is: 



S(q,w) 



S(q)M 



TTUJ [uj 2 p M /q 2 - M] 2 + [lUT] L ] 2 ' 



(2) 



where S(q) is the static structure factor. A Brillouin peak 
in the spectrum is thus expected close to 



lu la = q(M/p M ) 1/2 , 



(3) 



proportional to the longitudinal acoustic (LA) velocity 
cla = ^LA/q- The full width of the peak 



Tla = q 2 VL/pM 



(4) 



is proportional to the attenuation of the LA mode. In 
fact, the absorption of longitudinal acoustic modes is 
given by a = Y la/^-Cla- The integrated area of the spec- 
trum is proportional to the adiabatic compressibility xs- 
In the limit of low thermal conductivity here considered, 
thermal fluctuations contribute only to the low frequency 
region of the spectrum, giving rise to the Rayleigh peak 
centered at frequency zero and linewidth of about 0.1 
GHz. The total intensity, taking into account also this 
last contribution, is proportional to the isothermal com- 
pressibility xt- The ratio of the intensity of the Rayleigh 
line over the Brillouin lines (xt — Xs) /xs = 7 — 1 = Rlp 
is known as Landau-Placzek ratio. 

In molecular liquids the acoustic modes may couple 
both with the molecular internal degrees of freedom and 
with the a process, giving rise to additional mechanisms 
of acoustic damping. These relaxation effects are also re- 
sponsible for a peak at zero-frequency (Mountain peak) 
and can be taken into account by a generalization of the 
hydrodynamic equations, introducing an ^-dependent 



M{uj) = Moo + AM(w) + tutrix, 



(5) 



where Moo is the high-frequency (unrelaxed) longitudi- 
nal modulus and r/oo is the high-frequency longitudinal 
viscosity which accounts for relaxation processes char- 
acterized by a time so short to be out of the exper- 
imental frequency window. In the low frequency (re- 
laxed) limit, AM(w -> 0) = (M - Moo) + iuu{r] -ijoc), 
where Mq (rjo) is the zero-frequency longitudinal modu- 
lus (viscosity). In this limit, the complex longitudinal 
modulus becomes M(u —►()) = Mo + iuirjo. This is 
the previously described simple hydrodynamics regime, 
where 770 = Vl- I n the opposite limit, the unrelaxed one, 
AM(lu — » 00) =0 and the longitudinal modulus takes the 
simple form: M{uj) = + icur/oo- This coincides with 
the Voigt model of viscoelasticity which is used since long 
time to describe the acoustic damping in solids (see, for 
instance, |44j). In this limit the spectrum takes the same 
form of Eq. J5J , but with M replaced by M^ and t]l by 
rjoa . The intermediate condition between the relaxed and 
unrelaxed limits is the one where the relaxation processes 
strongly couple with the acoustic waves, and the detailed 
description of the shape of the Brillouin spectrum re- 
quires the knowledge of the frequency dependent part of 
the modulus, AM(u). In this condition Eq. (J2J can still 
be used to fit the spectrum in a narrow frequency region 
around the Brillouin peaks to obtain values for the "ap- 
parent" sound velocity and attenuation (see for instance, 
|2fl EEI]'!. i.e. their values at the frequency of the peak. 
This is the "traditional" way of analyzing Brillouin spec- 
tra, safe enough also for measurements obtained with low 
contrast and non-tandem interferometers. The analysis 
of the temperature behavior of these parameters gives 
useful information on the number and typology of relax- 
ation processes active in the Brillouin frequency window, 
and will be adopted in the first part of the data analysis, 
paragraph IIV Al 

In the most general case of coupling of relaxation pro- 
cesses with density fluctuations, the dynamic structure 
factor measured with both the BLS and the IXS tech- 
niques can be expressed in terms of the generalized mod- 
ulus of Eq. JSJ: 



S(q,u>) 



S{q)M 



M"(w) 



iTLU [LU 2 p M /q 2 - M'{lo)} 2 + [M"(uj)Y 



(6) 

In this case, the simplest possible situation is that of 
an exponential relaxation described by the Debye for- 
mula AM(u) = (Mo — Moo) /{I + iu>r), where r is the 
relaxation time. Notice that Mq < Moo, i-e. the modulus 
increases with increasing frequency. The imaginary part 
of the modulus shows a peak at w = t _1 corresponding 
to a maximum of acoustic energy absorption. It occurs, 
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for instance, when the frequency of the acoustic mode 
matches the rate of energy exchange with some internal 
degree of freedom. The spectral width of the peak of 
M"(u>) for a Debye relaxation is of 1.14 decades, in the 
same frequency region where the corresponding real part 
changes from Mq to M^ (dispersion). Viscoelastic sys- 
tems approaching the glass transition are characterized 
by a dispersion which extends in a much wider frequency 
range (2 or 3 decades is not a rare case). This stretch- 
ing phenomenon is frequently accompanied by a power 
law dependence of M"(u>) which can be described by the 
Cole-Davidson relaxation function 

AM(w) = (M ~ Moo)/(l + turf, (7) 

where (3 is the stretching parameter, being M"(ui) oc uj^ 13 
on the high frequency side of the peak in M"(u>). Fitting 
Brillouin spectra with this phenomenological ansatz for 
the longitudinal modulus gives the values of the ampli- 
tude (Mo — relaxation time and stretching param- 
eter at each measured temperature. The study of the 
temperature behavior of the a-relaxation parameters is 
the main topic of the physics of the glass transition. This 
is the content of Sec. Ivl 

No mention has been given up to now of the micro- 
scopic origin of the a-relaxation. MCT suggests that in 
the liquid above T gi the non-linear coupling of density 
fluctuations characterized by different values of q is the 
mechanism responsible for it. This theory gives a closed 
set of equations relating the dynamics of density fluctua- 
tions to the static structure factor of the same liquid. A 
full check of the predictions of MCT is beyond the aims 
of the present work. We limit ourselves to note that the 
Cole-Davidson relaxation function is usually adequate to 
represent the slow part of the relaxation function as ob- 
tained by the MCT [3 and its amplitude will be used 
in Section IV CI to check the existence of the square root 
singularity in the non-ergodicity factor, as predicted by 
the theory mEl. 



III. EXPERIMENTAL DETAILS 

A. Brillouin light scattering measurements 

Glycerol of 99.5% purity purchased from Aldrich 
Chemicals has been employed. The sample was fil- 
tered through a 0.22 fim teflon filter directly into the 
cylindrical pyrex cell (inner diameter w 1 cm) used for 
the measurements and slowly degassed under vacuum. 
The cell was sealed without breaking the vacuum. Bril- 
louin spectra were collected in a wide temperature range 
comprising the whole liquid - normal and supercooled - 
range (melting point, T m = 291 K) down to the glassy 
state (T g — 187 K). The sample was placed in a copper 
holder which was used to regulate the temperature. For 



the measurements above room temperature, the sample 
holder was heated with a NiCr wire using a conventional 
dc power supply. The temperature was controlled by a 
Thor Cryogenics 3010 II model temperature controller 
and measured by a PtlOO calibrated resistance placed 
close to the sample. The measurements were started af- 
ter the sample temperature was thoroughly equilibrated 
and the temperature fluctuations were kept within ±0.1 
K during the spectra recording. For the low-temperature 
measurements a Cryomech ST405 cryostat was used. The 
sample temperature, varied with a typical 0.5 degree/min 
rate, was measured with a 410-NN-l Silicon Diode Tem- 
perature Sensor. A Coherent-Innova 300 model Ar + laser 
was operated with a typical power of «300 mW on a sin- 
gle mode of the A Q =514.5 nm line. The light scattered by 
the sample was analyzed using a Sandercock-type (3+3)- 
pass tandem Fabry-Perot interferometer characterized by 
a finesse of about 100 and a contrast > 5-10 10 |46l l47j . 
The elastic scattering from a dilute water suspension of 
Latex particles (120 nm diameter) was used to deter- 
mine the instrumental resolution function R(lo). In the 
first part of our work, we aimed at the determination of 
the sound speed and attenuation of longitudinal acous- 
tic modes. For this purpose the peak position and the 
linewidth (FWHM) of the Brillouin components are the 
only relevant parameters, making it possible to focus on 
the frequency region of the Brillouin peaks, where the po- 
larized (Ivv) spectra replace in good approximation the 
isotropic (Iiso) ones. In this spirit, polarized Brillouin 
spectra were collected in the back-scattering (8 w 180°) 
geometry in a wide temperature range, namely 47-441 
K, using a narrow T-grid (about 5 K above room tem- 
perature, 10 K below room temperature). Different Free 
Spectral Range (FSR) values were used for the measure- 
ments above and below room temperature (20 GHz and 
10-^8.3 GHz, respectively). The integration time was of 
«10 s / channel and the effect of the dark counts (< 1 
count / s) was checked to be negligible. Typical polarized 
spectra covered the 0.8-^25 GHz frequency range. In a 
successive step we selected some temperatures above T g 
for a detailed analysis of the shape of the relaxation func- 
tion. We collected in the 210-400 K temperature range 
(namely, 400, 350, 333, 326, 295, 285, 270, 255, 240, 225, 
and 210 K) polarized (Ivv) an d depolarized (Ivh) spec- 
tra using two different FSRs of 13.6 GHz and 50.0 GHz 
in order to cover the 0.8-100 GHz frequency range. Each 
spectrum (polarized and depolarized) results from the 
overlap of the two spectra measured with different FSRs. 
The isotropic spectra Iiso have been then obtained by 
subtraction of the Ivv and Iyj j spectra following the 
procedure described in Ref. |45|. Typical Iiso spectra 
are presented in Fig. at some selected temperatures. 
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FIG. 1: Isotropic Brillouin light scattering spectra as a func- 
tion of temperature for normal and supercooled liquid glyc- 
erol. The lines through the experimental points are the best 
fits obtained using Eqs. 6 and 14. 
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FIG. 2: A typical IXS spectrum (bottom panel) acquired 
at T=350 K and q—2 nm _1 shown together with the BLS 
spectrum at the same temperature (top panel). The lines 
through the experimental points are the best fits obtained 
from the combined BLS-IXS analysis using Eqs. 6 and 14. 



analyzers. The total energy resolution is of w 1.5 meV. 
The IXS spectra are directly proportional to S(q, ui). IXS 
measurements have been performed in the 295-400 K 
temperature range, for q values in 2-11 nm _1 interval, 
with a g-resolution of about 0.4 nm^ 1 . Experimental 
data were normalized to the intensity of the incident 
beam and the contribution of the empty cell was sub- 
tracted from each spectrum. A typical IXS spectrum 
acquired at T=350 K and q=2 nm _1 is shown in Fig. [21 
together with the BLS spectrum at the same tempera- 
ture. 



C. Ultrasonic measurements 



B. Brillouin inelastic X-ray scattering 
measurements 

Inelastic X-ray measurements were carried out at the 
very high energy resolution IXS beamline (ID16) of the 
European Syncrotron Radiation Facility [48| using the 
Si (11 11 11) reflection for both the monocromator and the 



The temperature behavior of the relaxed sound veloc- 
ity Co was measured by an ultrasonic technique at 5 MHz 
[49j, i.e. at a frequency four orders of magnitude lower 
than that corresponding to the LA modes revealed by 
Brillouin light scattering. A PZT transducer was placed 
deep into the sample, at a distance L = 5mm from the 
bottom of a brass cell parallel to the surface of the trans- 
ducer. The system as a whole works as an acoustic Fabry- 
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Perot, so that the reflection coefficient, R, shows a series 
of minima corresponding to the condition L = to A/2, 
where A is the acoustic wavelength and m is an integer 
number. The reflection coefficient is measured as a func- 
tion of frequency by means of an HP8753A Network An- 
alyzer. The distance Av between two neighboring min- 
ima, the free spectral range, is thus given by the relation 
Av = c/2L, and the Fourier transform of R shows a well 
defined maximum at t = 2L/c, where c is the velocity of 
the acoustic wave in the sample. From the measurements 
of t, the velocity c of longitudinal acoustic waves was ob- 
tained at several temperatures in the range 306.5-389.5 
K. 



IV. DATA ANALYSIS 

A. Acoustic properties of glycerol 

The first step of our investigation is a traditional acous- 
tic analysis of the Brillouin spectra consisting in the mea- 
surement of the frequency position and linewidth of the 
Brillouin peaks as a function of temperature. It has been 
shown in Sec. Ill Al that the isotropic Brillouin spectrum, 
for frequencies close to those corresponding to the lon- 
gitudinal acoustic (LA) peaks, can be represented by a 
damped harmonic oscillator (DHO) model: 



I la (lo)=I° l 



LA h? 



LA 



+ [uT LA ]* 



(8) 



where lola/2tt and Tla/2tt correspond to the frequency 
position and to the full width at half- maximum (FWHM) 
of the LA peaks. 

The spectra collected in the 47-441 K temperature 
range were fitted using Eq. (|8J) convoluted with the ex- 
perimental resolution function. The values of ujla/2'k 
and Tla/^t: obtained with this procedure are reported in 
Table Q] The values of the (apparent) longitudinal sound 
velocity and of the (apparent) longitudinal kinematic vis- 
cosity have been derived using the relations: c = ujlaIi 
and D — Tla/q 2 , and are shown respectively in Fig. |3] 
and Fig. 01 where they are compared with literature data 
obtained from different techniques. Numerical values for 
c and D are reported as well in Table [I] The momentum 
transfer q at different temperatures has been calculated 
(see Tab. [J using the refractive index n and the rela- 
tion q = 2nki valid in the backscattering configuration. 
The T-dependence of n has been evaluated from that of 
the mass density p using the Clausius-Mossotti relation. 
Finally, the p values, in the temperature range here in- 
vestigated, have been obtained from the relationships: 



p{T > T g ) = 1.272 - 6.55 • 1CT 4 (T - 273.15), 



(9) 




p(T < T g ) = 1.332 - 3.20 • 10~ 4 (T - 187), (10) 



75 150 225 300 375 450 525 
T(K) 



FIG. 3: Sound velocity data from the present Brillouin light 
scattering experiment (full circles) compared with literature 
data: stars represent previous BLS measurements from |3^|. 
full squares represent ultrasonic measurements at 15 MHz 
from (23]. Limiting high frequency values obtained from our 
BLS full spectrum analysis are reported as open circles, to- 
gether with the expressions given in Rets, ib m m m 

(dash-dot line) . Limiting low frequency data determined from 
our US measurements are reported as open triangles tog ether 
with linear expressions given in Refs. [2'A . I'2 ll I2SL l3Cj (full 
line) . 



with T in K and p in g/cm 3 . Eq. © is from Ref. [24j] . and 
Eq. i|10fl has been obtained using the thermal expansion 
coefficient of the glass at 180 K, a g i ass = 2.4 • 10~ 4 K _1 

Looking at figures and 0] it can be seen that our 
data show the well-known features expected in presence 
of the structural relaxation process: a marked disper- 
sion of the sound velocity, accompanied by a maximum 
of the absorption. In particular, the hypersonic veloc- 
ity shows the typical S-shape dispersion curve (full cir- 
cles) that is bounded by the limiting high (unrelaxed) 
and low (relaxed) frequency values of the velocity, 
and Co, respectively The change in the slope of c(T) 
at about 187 K indicates the occurrence of the liquid- 
glass transition. We report in the same figure the val- 
ues obtained by ultrasonic measurements of the longitu- 
dinal modulus (squares) [2j| showing that a change of 
the probing frequency from the MHz to the GHz domain 
moves the dispersion region towards higher temperatures 
and broadens the temperature region where the a relax- 
ation is active. Fig. [3] also shows different linear extrap- 
olations of the unrelaxed adiabatic velocity taken from 

Refs. [UlllHil ( dash " dot line )> together with the val- 
ues of Coo obtained from the full spectrum analysis (open 



TABLE I: Peak position (lu la /2%) and FWHM (T la /2ty) of 
the Brillouin peak reported in the whole range of the investi- 
gated temperatures. These parameters have been determined 
by fitting the spectral region around the Brillouin peak us- 
ing the DHO model, Eq. JHJ, convoluted with the instrumen- 
tal resolution. The values of FWHM for temperatures lower 
than 196 K have not been reported because, on decreasing the 
temperature, the linewidths become too small, with respect 
to the resolution function, to be reliably extracted. Column 
2 reports the values of the exchanged momentum q. Columns 
4 and 6 report the longitudinal sound velocity and the longi- 
tudinal kinematic viscosity, respectively. 
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FIG. 4: Apparent kinematic viscosity Tla/i 2 derived from 
the analysis of the Brillouin spectrum (full circles) compared 
with literature data. Crosses are from Ref. [3(J, full triangles 
from Ref. [321 ]. Limiting high frequency values from Refs. ]3flL 
HUi l4lll are reported as diamonds, open triangles, and open 
circles respectively. Stars are the unrelaxed values determined 
from the BLS full spectrum analysis presented here. 



circles) which will be presented in the next section. Fi- 
nally, Fig. reports values of the relaxed sound velocity, 
whose temperature behavior is essential for the full spec- 
trum analysis of Sec. IIVBI In particular, several sets of 
data for Co are available in the literature |24l 1281 130| , ob- 
tained in different and narrow temperature ranges below 
360 K and above 440 K. They are typically described by 
almost linear trends with different slopes. In addition, we 
present here new ultrasonic data measured in an interme- 
diate temperature range between 306.5 K and 389.5 K. 
These data are magnified in Fig.[S]and the parameters for 
the linear approximation of cq in different temperature 
regions are reported in Tab. |n| Once these partial sets 
of data are shown all together in Fig. a globally non- 
linear behavior of cq(T) becomes clear, similar to that 
observed in many glass- formers |50|. Conversely, when 
we plot M = pcQ, as a function of T (inset of Fig. EJ, 
this quantity can be reasonably fitted by a straight line in 
a wide T-range from 245 to 505 K. The parameters of the 
fit are reported in Tab.HTl This linear behavior of Mq(T) 
is in agreement with what observed by Christensen and 
Olsen in a restricted temperature region [HH ]. 

Focusing now our attention on the longitudinal kine- 
matic viscosity, the well defined maximum close to 350 
K in Fig. 0] corresponds to the matching between the 
relaxation rate and the characteristic frequency of the 
longitudinal acoustic excitations. In the same figure we 
report as well the values of D determined by different 
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TABLE II: Temperature dependence of Co and Mo (T is in 
Kelvin) . 



Parameter 



Formula 



Units T — range Refs. 



c 2519-2.051 T 

c 2518-2.045 T 

c 2918-3.041 T 

co 2593-2.286 T 

Mo 8.110-1.183xl0~ 2 



T 



m/s 
m/s 
m/s 
m/s 
GPa 



232-343 K 
200-360 K 
440-534 K 
307-390 K 
245-505 K 



f24] 
[28] 
f30l 



"Present work. 



techniques. Particularly relevant is the comparison with 
x-ray scattering data. Generally, the IXS probing fre- 
quency is so high that, especially for low temperatures, 
the structural relaxation is definitively out of the spectral 
window. Therefore, IXS typically probes the unrelaxed 
values of T/q 2 , i.e. r/^/p. Fig. 0] shows that the kine- 
matic viscosity data obtained by BLS and IXS are con- 
sistent, within their error-bars, in a range of about 80 K 
below T g beside the fact that they correspond to q 2 val- 
ues which differ by a factor 10 3 . Below T g , the structural 
relaxation is located at frequencies lower than one Hz, 
and we cannot expect any significant contribution from 
it to the acoustic loss in the GHz region. The absence of 
an excess damping in BLS data over the IXS unrelaxed 
ones at T g points in favor of an a relaxation-only scenario 
for the dynamics of glycerol. This is a peculiar behavior 
of glycerol, and, possibly, of associated liquids, different 
from that observed in several fragile systems like OTP, 
DGEBA and PB H3, where the presence of an excess 
attenuation in the GHz range reflects the coupling of the 
LA modes with internal, molecular degrees of freedom, 
giving rise to fast relaxation processes of thermal nature. 
Below 100 K, the BLS values of D [^3 start to decrease, 
suggesting a scaling of T with a power of q higher than 2. 
A dynamic contribution to the sound attenuation hav- 
ing influence on low frequency density fluctuations has 
been recently proposed as the origin of this phenomenon 

mil. 

Additional information on the number and nature of 
the relaxation processes active in glycerol can be obtained 
from the high temperature-low frequency non dispersive 
region. This regime is well described by simple hydrody- 
namics |l4j , and the acoustic absorption, in low thermally 
conducting liquids, is given by: 



OJ 



VL 

2pc 3 



(11) 



where rj L = rj b + 4/3^s and tjl, ?y b , and r/s are the longi- 
tudinal, bulk and shear viscosity, respectively. In simple 
liquids, where rjb is small, rjs gives the most important 
contribution to the acoustic attenuation, giving rise to 
the so-called classical absorption [Til ]: 
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FIG. 5: Relaxed sound velocity. Our Co measurements col- 
lected at 5 MHz are reported as full circles and are com- 
pared with previous determinations: triangles refer to mea- 
surements reported in Ref. [24lj2p and stars represent the 
linear expression given in Ref. |3fl|. This whole set of data 
does not show a linear behavior. However the same data, 
once reported in terms of relaxed longitudinal modulus, Mo, 
do show a linear behavior, as reported in the inset of the fig- 
ure. The linear representation of Afo is reported in terms of 
Co as dotted line. 



From Eq. i|ll|) and Eq. i|12fl , the relationship between the 
longitudinal and the shear viscosity can be easily written 
as: 



'Is 



4 a 
3 a c i 



(13) 



A different temperature behavior of the a J a r i ratio 
distinguishes Kneser and associated liquids [lj]. The 
Kneser-type systems, as for example ionic liquids and 
the majority of organic liquids, show a temperature de- 
pendent ajoLci ratio ranging between 3 and 400. In these 
liquids, the excess absorption over the classical value is 
due to exchange of energy between internal (rotations 
and vibrations) and external (acoustic modes) degrees 
of freedom, showing a marked temperature dependence 
[Til ]. On the contrary, for the associated liquids, often 
characterized by CH3, CH2, and OH groups, the bulk 
viscosity is essentially only related to the structural re- 
laxation, the contributions coming from the internal de- 
grees of freedom being negligible. In these systems, the 
acoustic measurements performed in the non-dispersive 
region show a temperature independent ratio a/a c i lying 
between one and three [T3 |. In particular, for glycerol, 
it has been proven that the aja c i ratio is approximately 
constant and equal to 1.78 |14| . 
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FIG. 6: Shear viscosity determined from our Brillouin data 
as r]s = r\i,l1 A = Y LAp/IAq 2 (full circles). These values 
are compared with those obtained from direct shear viscosity 
measurements: triangles refer to measurements reported in 
Ref. |H^|. stars represent measurements from Ref. |55|. and 
the straight line indicates the expression reported in Ref. (24) . 



On this basis, starting from the values of ijl obtained 
from the present BLS data, and assuming the same tem- 
perature dependence for shear and volume viscosity, rjs 
has been obtained using the Eq. (T3J), and a/a c i = 1.78 
[l4| . In Fig. we report such an evaluation for ns to- 
gether with viscosimetric measurements pil 155) . In 
the high temperature region, the two data sets show a 
common trend, reinforcing the hypothesis that no fast 
thermal relaxations are responsible of the dynamics of 
glycerol, at least in the frequency and temperature re- 
gions covered by Brillouin light scattering. This state- 
ment is even more clear if Fig. [S] is compared with Fig. 
7 of Ref. H3, where D(T) for the Kneser liquid OTP 
is reported. In this latter case, a marked temperature 
dependent a/a c i ratio has been observed, attributed to 
the effects of a thermal relaxation, whose presence has 
been recently supported by molecular dynamics simula- 
tions H2- 



B. Full spectrum analysis 

The traditional analysis, presented in previous section, 
helps us to choose an ansatz for the frequency dependent 
longitudinal acoustic modulus M*(u>). To this respect, 
the most relevant indications we gained from this ap- 
proach are: i) at temperatures close to the glass tran- 
sition the (apparent) kinematic viscosity obtained from 
BLS is consistent with the corresponding unrelaxed val- 



ues determined by IXS, and it) in the high temperature 
limit the shear viscosity estimated from BLS collapses on 
that obtained by viscosimetric measurements. Both these 
results indicate the absence of thermal relaxations in 
glycerol, in agreement with what inferred some 40 years 
ago on the basis of ultrasonic measurements . There- 
fore, in order to model M*(u), we may suppose that 
the structural relaxation is enough to catch the leading 
contributions to the relaxational dynamics of glycerol in 
the GHz frequency window. Relaxation processes so fast 
to be out of the experimental frequency window are de- 
scribed by an instantaneous term in the modulus function 
and their effects are taken into account by the unrelaxed 
viscosity r/^,. Following these indications, our ansatz for 
the M*(lo) function consists of a Cole-Davidson function 
for the a relaxation plus an instantaneous contribution: 



M (oj) = Moo + -jz — ; — : 773- + «w?7o 

(1 + lWT a )P 



(14) 



Substituting this expression for M*(u>) into Eq. 6, we 
have obtained a model function to fit to the BLS spectra. 
In particular, this full spectrum analysis has been carried 
out for temperatures higher than T g , adopting the follow- 
ing procedure: (i) The thermal diffusion contribution is 
considered as a Dirac 5 term centered at zero frequency, 
whose amplitude is proportional to that of the adiabatic 
part of the spectrum via the LP-ratio. In fact the charac- 
teristic time of the thermal diffusion mode as measured, 
for instance, in Ref. [3l| and properly rescaled to the q 
values accessed to by BLS is located at about 10 ns, i.e. 
more than three decades far from the Brillouin peak in 
the whole temperature region considered here. Moreover, 
even at the highest investigated temperatures, if the ther- 
mal contribution is explicitly considered as a Lorentzian 
peak centered at zero frequency, the quality of the fit 
is unchanged together with the values of the fitting pa- 
rameters, (ii) The values of the relaxed adiabatic sound 
velocity Co have been fixed to those reported in Tab. II I II 
(iii) r)oo/ p has been left to vary in the fit ha ving as start- 
ing point the value determined by IXS p39l l40il4ll| . The 
weak T-dependence and the experimental uncertainty of 
IXS data discouraged us from fixing r)oal P to a constant 
value, (iv) The unrelaxed adiabatic sound velocity Coo 
has also been taken as an adjustable parameter in the 
whole T-range, while T a , and (3 have been considered 
free or fixed depending on the temperature region inves- 
tigated. In particular, in the low temperature region (T 
< 255 K), being the structural relaxation out of the Bril- 
louin window, r a was fixed and (3 left free since it affects 
the shape of the Mountain region of the spectrum. In 
this region t q has been imposed to be proportional to 
the shear viscosity 77s, i.e. r Q = J rjs, where J is a con- 
stant coefficient havin g th e dimensions of a compliance 
(see for example Ref. |2(|). In contrast, in the highest 
temperature region (T > 270 K) the structural relax- 
ation becomes progressively active in the GHz frequency 
window, so that we are able to obtain r a from the fit, but 
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we need to fix the value of (3. Its value has been chosen 
within the limits suggested from Ref. |2^|, i.e. 0.40±0.05, 
and also consistent with our results at lower T. For some 
selected temperatures, namely 295 K, 350 K, and 400 
K, we simultaneously fit BLS (g=0.036 nm _1 ) and IXS 
(q—2 nm -1 ) spectra, following the combined fitting pro- 
cedure described in Ref. |18| . In this procedure we use 
a single IXS spectrum collected at q=2 run -1 , since the 
spectra at higher g's progressively depart from the hydro- 
dynamic regime. In the joint-analysis, the values of 
and 7700 are mainly related to the position and linewidth 
of the Brillouin peaks of the IXS spectra, while r a and j3 
are mainly related to the shape of BLS spectra. The ob- 
tained best fitting curves are shown in Fig.^and Fig.Elas 
full lines. The whole set of fitting parameters is reported 
in Tab. IIIII The relaxation times obtained from the BLS 
data-analysis are also shown in Fig. together with sev- 
eral literature data. In particular, in the same figure 
we report three different evaluations of r Q for T=350 K: 
the value obtained from the described fit procedure (full 
circle); the value of r max (down triangle) corresponding 
to the maximum of the Cole-Davidson function used to 
model the BLS spectra; and the characteristic time of 
longitudinal acoustic mode tla = ^-/^vla (star), i.e. 
the time corresponding to the maximum of the BLS ab- 
sorption data in Fig.0] We underline that this last eval- 
uation of r Q is model independent. Fig. [7| shows that, 
at T=350 K, r max and tla coincide. This result repre- 
sents a strong support of the reliability of the values of 
the structural relaxation time as deduced from the full 
spectrum analysis presented here. 



V. DISCUSSION 

A. The structural relaxation time and the 
stretching parameter 

The temperature behavior of the structural relax- 
ation time determined in the previous paragraph is 
compared with that obtained by different spectroscopic 
techniques (dielectric spectroscopy, impulsive stimulated 
Brillouin scattering, impulsive stimulated thermal scat- 
tering, ultrasonic techniques, depolarized light scatter- 
ing) in Fig. The inset of the figure shows that the 
data obtained by these techniques, once rescaled by a 
constant, collapse into a single master curve, indicating 
that the a-scale universality is obeyed in glycerol in a 
wide time-temperature range. A simple explanation for 
the difference in the absolute values of r is the different 
microscopic correlators measured by the different tech- 
niques. Next in importance is the definition of the char- 
acteristic time itself. In some cases, like in BLS and 
ultrasonic measurements, the longitudinal acoustic mod- 
ulus formalism has been used, with a specific CD complex 
function to describe the structural relaxation process. 
In some other case, like in dielectric spectroscopy mea- 
surements, the compliance formalism was chosen and, al- 
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FIG. 7: Structural relaxation time of glycerol from the BLS 
full spectrum analysis compared with various literature data. 
Full circles and squares represent our BLS values obtained in 
the fitting procedure by leaving the r a parameter free or by 
imposing it to be proportional to the shear viscosity, respec- 
tively. The star represents the value tla corresponding to 
the maximum of the acoustic attenuation for the longitudinal 
modes; the down triangle marks the value of r a at 350 K cor- 
responding to the maximum of the Cole-Davidson function 
used to model the BLS spectra. The open squares represent 
ultrasonic data taken at 2, 10, and 27 MHz |25|. pluses are 
from Ref. |3lj ] , crosses are from Ref. [5(3 ] , open circles are from 
Ref. 1371. The short-dashed line is the polynomial law given 
in Ref. |28q . The dash-dot line is the Vogel-Fulcher equation 
reported in Ref. I2H describing photon correlation data. The 
solid line is the Vogel-Fulcher equation obtained in Ref. |2fi| | 
from a fit to ISTS data. In the inset, the same data are re- 
ported after rescaling by a constant. The collapse of all the 
data on a single master-curve indicates the validity of the 
a-scale universality in glycerol over a wide time-temperature 
region. 



though a CD relaxation function was used, its parameters 
are not analytically related to those of the CD modulus 
function. Moreover, some authors reports the values of 
T a directly obtained from the fit p^ . l28ll5r3 |: others report 
the average value < t > (/3r, in case of a CD function; 
(r//3)r(l//3), where T is the Gamma function, in the case 
of a Kolhraush- Williams- Watt function) [2^, |^| ; others 
report the value corresponding to the maximum of the 
absorption T max |25l l37| . 

Despite these technical obstacles, one reasonably ex- 
pects that, at least in those cases where the same corre- 
lator is measured, a quantitative agreement between re- 
laxation times obtained by different techniques should be 
reached. In particular BLS, ultrasound, ISBS, and ISTS 
all probe density fluctuations. The agreement between 
BLS, ultrasonic, and ISBS r-modulus data (present work, 
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TABLE III: Parameters of the full spectrum analysis of liquid and supercooled glycerol at the temperatures reported in column 
1. Columns 2, 3, and 4 report respectively the characteristic time of the a-CD function, the stretching parameter, and the 
average value of r determined as (r a ) = /3r Q . Columns 5 and 6 report the high and low limiting values of the sound velocity, 
respectively. Column 7 reports the instantaneous kinematic viscosity contribution. 
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"Fixed on the basis of the temperature behavior of the shear vis- 
cosity. 

'Free parameter in the fit. 

c Ref. (24) . 

d Present work. 

c Value obtained from the interpolation of our data of co and those 
reported in Ref. |30l . 



IU, [H) is evident in Fig. □ Different from this, ISTS 
r-compliance data |29l l31| are shifted by about a factor 
ten towards longer times. Even after a conversion of the 
BLS data to the compliance formalism, a shift of about a 
factor three persists between the BLS and the ISTS data. 
At present, we cannot find any explanation for this in- 
teresting anomaly, and we limit ourselves to suggest to 
perform this comparison on other glass-forming systems. 

The present full spectrum analysis provides informa- 
tion also on the T-dependence of the shape parameter 
Pcd of the a-relaxation function. In particular, an al- 
most constant value, namely 0.37±0.03, was found for 
temperatures between 270 and 400 K, close to the value 
0.40±0.05 obtained in a previous ultrasonic investiga- 
tion |25| . The value of the stretching parameter de- 
creases with temperature, as shown in Tab. IIIII - the 
same trend being also found in the dielectric investi- 
gation reported in Ref. the only one which cov- 
ers an extremely wide temperature-frequency range. For 
a quantitative comparison with the values obtained by 
other techniques where compliance relaxation functions 
are used, like ISTS, photon correlation spectroscopy, and 
specific heat spectroscopy, we converted our modulus- 
based data into compliance-based data obtaining the 
values reported in Tab. IIVI Moreover, for those tech- 
niques where the stretching parameter was referred to 
a Kohlrausch- Williams- Watt function (KWW), we have 
transformed 0kww into @cd using the equations re- 
ported in Ref. 57]. Our BLS data compare fairly well 
with an average value of f3cD~0-^5±0.05 estimated from 
ISTS measurements at different T and exchanged mo- 



TABLE IV: Summary of the shapes of different compliance 
functions. 



Technique 


PCD 


Temperature range 


Refs. 


BLS 


0.40±0.05 


210 K 


present work 


BLS 


0.40±0.05 


225 K 


present work 


BLS 


0.53±0.05 


240 K 


present work 


BLS 


0.56±0.05 


255 K 


present work 


ISTS 


0.45±0.05 


200-340 K 




PCS 


0.40±0.03 


195-225 K 


[26, 27] 


Ultrasonics 


0.42±0.05 


235-300 K 


[25J 


Specific heat 


0.51±0.03 


195-225 K 


f58] 


Dielectrics 


0.55-^0.65 


195-300 K 


f37] 



mentum, /3 CjD =0.40±0.05 from PCS data jH H3 and 
/3<2p=0.51±0.03 from specific heat spectroscopy data 



B. The limiting high frequency sound velocity 

The limiting high frequency values of the adiabatic 
sound velocity obtained from the fit of the Brillouin spec- 
tra are shown in Fig. |H1 In the same figure, we report 
two different ansatzs for the unrelaxed speed of sound 
obtained assuming either a linear T-dependence of My 
or a linear T-dependence of c x %MM%. In Ref. [loj|, 
using stimulated Brillouin gain spectroscopy, a slight but 
continuous curvature in the data was observed and 
two different expressions were reported. Our values of 
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Coo agree with those given in Ref. [3(j for the 185.5-281.4 
K temperature range (stars in Fig. • However, we ob- 
serve that an extrapolation of this equation to higher 
temperatures is not able to reproduce our BLS results. 
In the cases of Refs. 0> l^sj - the linear expressions re- 
ported in Fig.|Hlare derived as extrapolations of the data 
collected at low frequencies, in a region where, possibly, 
complete unrelaxed conditions are not fulfilled. The val- 
ues of Coo obtained from the present BLS analysis display 
a slope smaller than in previous determinations and a de- 
creasing trend on increasing temperature with a possible 
change of regime at about 310 K. This estimation is in 
remarkable agreement with a previous IXS determina- 
tion, namely T x = 300 ± 20 K [4£| , which was suggested 
to mark the microscopic transition between two different 
dynamical regimes related to the glass and liquid phases, 
respectively. 

It is interesting to notice that a different representa- 
tion can be given for the same data by plotting the lon- 
gitudinal modulus Moo = pCoo as a function of T. The 
inset of Fig. [S] shows that Moo(T) can be reasonably ap- 
proximated by the linear behavior Moo(T) = [26.2(2) — 
0.05(1)T] GPa. This is similar to what found for the 
T-dependence of the limiting shear modulus of glycerol 
[59j . This linear behavior of the unrelaxed moduli is not 
universal in supercooled liquids. Indeed, in many glass- 
formers, a linear extrapolation leads to unphysical low 
values of Moo m the relaxation region. Rather, the limit- 
ing compliance Joo , the reciprocal of Moo > is found often 
to vary linearly with temperature [2J, |5pJ . 



C. The non-ergodicity factor 

A detailed analysis in terms of the predictions of the 
mode-coupling theory is beyond the subject of this work, 
we only present here an estimation of the non ergodicity 
parameter, /, derived from our isotropic Brillouin spec- 
tra. With BLS, in fact, we directly probe the density- 
density correlators, and thus we can extend the tem- 
perature interval recently investigated in Ref. [2^|. In 
the zero exchanged momentum limit, q — > 0, the non- 
ergodicity factor, fq, can be calculated via the equation 
fo = 1 — Cq/c^ jijj. This is equivalent to integrate the 
Mountain peak in quasi elastic light scattering spectra 
or the quasi elastic peak in neutron scattering spectra. 
Therefore, the fo parameter has been calculated from 
the Co and Coo values reported in Tab. 11111 The temper- 
ature interval investigated here, namely 210-400 K, in- 
cludes several previously different T c estimations. Look- 
ing at Fig. |21 it is evident that fo decreases with increas- 
ing temperature, but no signature of a cusp-like behavior 
emerges in glycerol. Fig.[|J]also shows that our results are 
consistent with those recently obtained by Paolucci and 
Nelson [2{| who measured the non-ergodicity factor from 
the strength of the structural relaxation using the ISTS 
technique in the 228-268 K temperature range. However, 
although the existence of the critical temperature cannot 
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FIG. 8: Unrelaxed sound velocity. The values of Coo deter- 
mined from the BLS full spectrum analysis (full circles) are 
compared with various proposed extrapolations: dashed, dot- 
ted, and dash-dot lines are from Refs. |28|. |24|. and |3(J, 
respectively. Stars represent data of Coo from the linear ex- 
pression given in Ref. 13 ( 1 in the 185.5-218.4 K temperature 
range. Solid lines are a guide for the eye. In the inset circles 
represent the unrelaxed longitudinal modulus M x = pc^. 
The solid line represents the equation obtained by a linear fit 
of the modulus data: Moo(T) = [26.2±0.2-0.05±0.01 T] GPa. 

be supported here, we consider this check only a prelimi- 
nary test of the MCT predictions, and it is not possible to 
exclude the validity of a "non-idealized" mode-coupling 
description for glycerol. 



VI. CONCLUSIONS 

Brillouin light scattering measurements of glycerol 
have been presented, performed in a wide temperature 
region ranging from the high temperature liquid to the 
glassy state. Our investigation was aimed at understand- 
ing the dynamics of density fluctuations in a prototyp- 
ical intermediate glass-former in the high frequency re- 
gion where the precursor of the glass transition, i.e. the 
structural relaxation, is present above T g . The associated 
nature of glycerol plays in favor of a simple relaxation 
pattern, where the additional effects of internal thermal 
relaxations is negligible in the GHz region. Focusing on 
the narrow frequency region of the Brillouin peak, the 
temperature behavior of the characteristic frequency and 
lifetime of the longitudinal acoustic modes has been an- 
alyzed. In particular, we have drawn the attention on 
the behavior of the apparent kinematic viscosity that, at 
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FIG. 9: The non-ergodicity factor in glycerol as a function 
of temperature. The values of fo determined from our BLS 
full spectrum analysis are compared with the data collected 
at different wavevectors and reported in Ref. [29|. 



temperatures close to and below T g , approaches the un- 
relaxed value estimated from IXS data. Different from 
glass forming systems that show secondary relaxations 
of intramolecular nature, no excess of attenuation for the 
longitudinal acoustic modes has been revealed near the 
glass transition temperature. Moreover, the longitudinal 
to shear viscosity ratio is found to be almost tempera- 
ture independent in the non-dispersive high-temperature 
region. These findings suggest that an a-relaxation only 
model plus an instantaneous viscosity term are able to 



catch the leading contributions to the dynamics of the 
density fluctuations. The choice of this model was also 
supported by a combined analysis of light and inelastic 
x-ray scattering spectra. Starting from this guess, we 
analyzed high-resolution BLS spectra over two decades 
in frequency (full spectrum analysis) and we were able 
to obtain details on the relaxation parameters from the 
shape of the Brillouin peaks and of the Mountain region 
of the spectrum. 

The temperature behavior of the structural relaxation 
time, T a (T), obtained by this fitting procedure, has been 
compared with the results coming from different spectro- 
scopic techniques. The a-scale universality is obeyed for 
the structural process, i.e. the values of r a are propor- 
tional to those revealed by different experimental tech- 
niques. Moreover, a significant difference of the absolute 
values of our r a data is observed with respect to those 
obtained by ISTS, a technique which probes the density 
fluctuations as well. At the moment, we have no ex- 
planation for this anomaly. It would be interesting to 
investigate whether this anomaly holds for other glass- 
formers as well. Finally, preliminary considerations have 
been presented on the predictions of the mode coupling 
theory. In particular, in the analysis presented here, the 
calculated non-ergodicity factor, f , is found to decrease 
with increasing temperature, but no signature of a cusp- 
like behavior is revealed in glycerol. We can thus argue 
that the predictions of the idealized MCT are not fulfilled 
in this system. Of course, this does not imply that MCT 
is not appropriate to describe the dynamics of glycerol. 
Instead, this suggests that a more complete version of the 
theory should be applied in this case. To this respect, the 
model introduced by Franosch et al. 60] , including both 
linear and quadratic coupling contributions and tested on 
depolarized light scattering spectra of glycerol, seems to 
be a suitable candidate to be tested on spectra of density 
fluctuations. Work along this line is in progress. 
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